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Summary. - Previous work in the literature has studied the Hamiltonian structure of 
an R 2 model of gravity with torsion in a closed Friedmann-Robertson- Walker universe. 
Within the framework of Dirac's theory, torsion is found to lead to a second-class primary 
constraint linear in the momenta and a second-class secondary constraint quadratic in the 
momenta. 

This paper studies in detail the same problem at a Lagrangian level, i.e. working on 
the tangent bundle rather than on phase space. The corresponding analysis is motivated by 
a more general program, aiming to obtain a manifestly covariant, multisymplectic frame- 
work for the analysis of relativistic theories of gravitation regarded as constrained systems. 
After an application of the Gotay-Nester Lagrangian analysis, the paper deals with the 
generalized method, which has the advantage of being applicable to any system of differ- 
ential equations in implicit form. Multiplication of the second-order Lagrange equations 
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by a vector with zero eigenvalue for the Hessian matrix yields the so-called first-generation 
constraints. 

Remarkably, in the cosmological model here considered, if Lagrange equations are 
studied using second-order formalism a second-generation constraint is found which is 
absent in first-order formalism. This happens since first- and second-order formalisms are 
inequivalent. There are, however, no a priori reasons for arguing that one of the two is 
incorrect. First- and second- generation constraints are used to derive physical predictions 
for the cosmological model. 



PACS numbers: 04.20.Cv, 04.50, 04.60.+n, 98.80.Dr 
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1. - Introduction. 

In recent years part of the theoretical-physics community has started appreciating 
that constrained systems may be also studied using a tangent-bundle formalism [1] , rather 
than working on phase space. From the point of view of relativistic theories of gravitation, 
this approach appears relevant since it may lead to a fully covariant treatment of the 
gravitational field regarded as a constrained system, without having to use a 3+1 split of 
the Lorentzian space-time manifold. 

Relying on previous work on cosmology in [2], we perform the analysis, in a closed 
Friedmann-Robertson- Walker (hereafter referred to as FRW) universe, of a model of grav- 
ity with non-vanishing torsion where the Lagrangian density is proportional to the square 
of the scalar curvature. Of course, this is only a toy model, but its analysis deserves 
careful consideration because if one tries to quantize general relativity (GR) or theories 
with non-vanishing torsion within the framework of perturbative renormalization, one finds 
that quantized GR is perturbatively non-renormalizable, and in both theories the effective 
action acquires terms quadratic in the full Riemann-curvature tensor, i.e. terms propor- 
tional to R a b c dR abcd i R a bR ab , (-^°a) • I n the most general case, these theories, which are 
non-linear in the curvature, are either non-unitary or perturbatively non-renormalizable 
[3]. Moreover, they are much studied as gauge theories of gravitation [4-6], although the 
corresponding Cauchy problem at the classical level may be ill-posed. However, to improve 
our understanding of all difficulties and limits of the perturbative-renormalization program 
in quantum gravity, and to check the consistency or inconsistency of inflationary-universe 
models based on these non-linear theories, it appears essential to complete a detailed anal- 
ysis of these models of gravity. 

In the torsion-free relevant cosmological application was studied, for example, 

in [7], where an action functional was considered which is purely quadratic in the trace of 
the Ricci tensor. The motivation for that quantum-cosmological analysis was to obtain a 
Wick-rotated path integral whose integrand does not blow up exponentially after suitable 
confer mal rescalings of the 4- metric. As a first step of the program motivated so far, 
a similar model was proposed in [2] (see corrections in [8]), from a Hamiltonian point 
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of view, by applying Dirac's theory of constrained Hamiltonian systems at the classical 
level. The resulting idealized cosmological model is a constrained system with second- 
class constraints, which arise, as shown below, by virtue of torsion. In the following part 
of this introductory paragraph, to help the reader, we summarize the Hamiltonian analysis 
appearing in [2]. 

In our FRW model, spatial homogeneity and isotropy imply that, in a coordinate 
frame, the only non-vanishing components of the torsion tensor are S X q = S 20 2 = S 30 s = 
Q(t). Denoting by N the lapse function, by a(t) the 3-sphere radius of the closed FRW 
metric, and defining a = log(a), ^ = y, r = J N dt, ' = J^, the action functional 



and \x is a proportionality constant, denoted by /i° in equation (2.13) of [2]. Note that the 
point Lagrangian (1.2) is obtained by inserting the FRW hypothesis into the field-theory 
Lagrangian of (1.1). By a direct calculation, one can check that this procedure is correct in 
the particular case of FRW cosmologies in general relativity [9]. For alternative theories, 
such a method was found to be correct in [7] in the FRW torsion-free case. We are thus 
using the assumptions about the FRW symmetry to reduce ourselves to the study of point 
Lagrangians without having to build the complete Hamiltonian treatment for a generic 
curvature-squared field theory with torsion. 

Since the addition to the Lagrangian of a total derivative leads to an equivalent set of 
field equations, we use this property to cast the theory in Hamiltonian form. To eliminate 
the square of second-order derivatives appearing in L one is thus led to define [2,7,10] 




(1.1) 



after integration on the 3-sphere, takes the form I = J L dr, where 



L = /je 3a e~ 2a + 2(a') 2 - Ylya! + 16y 2 + a" - 2y' 



(1.2) 




(a' 



2y)z 



(1.3) 
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where (a' — 2y) is proportional to the trace of the extrinsic-curvature tensor of the 3- 
sphere, and z is obtained differentiating L with respect to the highest derivative, i.e. 

defining z = J^, where x = a — 2 J y dr. Thus, setting y = u', the Lagrangian defined in 
equation (1.3) becomes (cf. equation (2.1)) 

z 2 

L = 16z(u') 2 + 2z{a') 2 - 12 zu' a' + 2u'z' - z'a' + ze~ 2a - — e~ 3a . (1.4) 

For comments on this choice of variables, see again [2,7]. Hence, defining p a = J^-, 
p u = p z = -|p- (cf. section 2), one finds the primary constraint rs 0, where the 
weak-equality symbol ~ denotes an equality which only holds on the constraint manifold 
[11-13], and 

</>i = 2p a + p u - Azp z . (1.5) 

The corresponding effective Hamiltonian H on the whole phase space is given by H = 
H c + 70i, where H c , the Legendre transform of L, takes the form 

9 — In 

Hc = -4zpl + ^ + Z -^--ze-^ . (1.6) 

The constraint 4>i is preserved in time by requiring that its Poisson bracket with H, as 
defined in [11], should vanish: ps 0. This leads to the secondary constraint 

<p 2 = I6zp 2 z - 2p uVz + ^e" 3a - 8ze~ 2a . (1.7) 

The constraints (pi and 4> 2 are second-class, since, after evaluating their Poisson bracket, 
one obtains a function which does not vanish when (pi and (f> 2 are set to zero. Hence 7 



(t>2, H c 

can be obtained as 7 = —4 y as in equations (2.24)-(2.25) of [2]. Note that, since (pi 

and (p 2 are second-class, we can define a new Poisson bracket, the so-called Dirac brackets, 
in which second-class constraints can be set strongly to zero [11-13], i.e. they behave as 
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Casimir functions for Dirac brackets. This implies that the canonical Hamiltonian H c in 
equation (1.6) may be also written as 



H c = -2zp\ - p a p z + 



— ze 



-2a 



(1.8) 



4^ 



which formally coincides with the torsion- free result [7]. The effective Hamiltonian H, 
however, is not a linear combination of constraints and hence does not vanish in general. 
One then finds the field equations (2.27)-(2.32) of [2]. 

After this introductory paragraph, we can summarize the plan of our paper as follows. 
Section 2 performs a Lagrangian analysis of our second-class constrained system. The 
kernel of the pre-symplectic two-form, its non-vertical part, the secondary constraint $2, 
the second-order vector field solving the field equations and tangent to the constraint 
manifold are derived in detail. Section 3 presents instead a constraint analysis within the 
framework of the recently proposed generalized method. Results and open problems are 
presented in section 4. 

2. - Gotay-Nester Lagrangian analysis. 

The Lagrangian L of the model outlined in section 1 is more conveniently re-written 
for our purposes in the form 



where a', u', z' have been replaced by the tangent-bundle fibre coordinates v a , v u , v z 
respectively. The corresponding GN Lagrangian analysis [14] is as follows [15]. We first 
evaluate the Cartan one- form 6l = §^r dq l , and the pre-symplectic two-form ujl = dd^ 
[16]. Given a vector field Y belonging to the tangent bundle, and setting to zero the con- 
traction iy^L, one thus finds the kernel ker ujl of the pre-symplectic two- form. Moreover, 
denoting by El the energy function, with corresponding one- form dE^, the vanishing of 
the contraction i~ dE^ defines the secondary constraint $2 5 for Yq belonging to the non- 
vertical part of ker ujl- Denoting by Y an element of ker u>l, this Lagrangian definition 



L = I6zv^ + 2zv^ — 12zv u v a + 2v u v z — v z v a + ze 



-2a 




(2.1) 
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of constraints is clearly understood acting with i Y on both sides of the Euler-Lagrange 
equations written in the form 

i r uJ L + dE L = , 

and using the identity iyiv^L = — iv^Y^L- This yields the condition = —iydE^. Thus, 
unless iydEL is identically vanishing, such calculation shows that $2 = iy c dE^ is actually 
the secondary constraint of the theory (section 4). This constraint is then preserved by 
requiring that its Lie derivative along the vector field V which solves the Lagrange field 
equations should vanish. 

Indeed, from equation (2.1) one easily finds that the Cartan one-form and the pre- 
symplectic two-form are given by 

9l = (&zv a — 12zv u — v z ^jda+^S2zv u — 12zv a + 2v z ^jdu+(2v u — v^dz , (2.2) 



(Av a — 12v u ^jdz + 4z dv a — 12z dv u — dv z ~^ A da 
(^32v u — 12va^jdz — 12z dv a + 32z dv u + 2dv z ~^ A du 
2dv u - dva\ A dz . (2.3) 



u L = 
+ 
+ 



Moreover, for a vector field Y of the tangent bundle, whose general decomposition is 
Y = Y a ( q ,v)± + Y.{ q ,v)± + Y.( q ,v)± 

the contraction with ujl is evaluated according to the formula 



7 



Lagrangian theory of constrained systems: cosmological application 
After a lengthy calculation, this yields 



i y ool = 

+ 
+ 



4zY Va - \2zY„ - Y v + ( 4v a - 12v u )Y 2 



(Ave, - I2v^j: 



da 



- 12zY Va + 32zY Vu + 2Y Vz + (s2v u - 12v a )Y 2 



du 



- Y Va + 2Y Vu + (l2v u - 4 Vc ^Y a +(l2v a - S2v u ^Y u jdz 



+ 
+ 



AzY a + \2zY u + F 2 ] dv a + 



12zY a - 32zY u - 2Y Z dv u 



(Y a - 2Y u ^j 



dv 2 



(2.6) 



Thus, if we set to zero iy^L-, the vector field Y G ker ujl is found to take the form 



Y = Y a *- + 1 -Y a *- 
da 2 du 



d d 
2zY a — + Y Vc , —+ 



dz 



1 



Y Va + (2v u -v a )Y c 



_d_ 



+ 



- 2zY Va + Azv a Y a 



_d_ 

dv 2 



(2.7) 



Note that, in equation (2.7), Y a and Y Va remain arbitrary functions of a, u, z, v a , v u and 
v z . For Y a = we get the vertical kernel. For them dE^Y) = is identically satisfied. By 
setting Y Va = we get vector fields Yc giving rise to constraints. Thus, since the energy 

dL 



function El = v l M^ — L is, in our case, such that 



E L = 16zv u + 2zv a - 12zv u v a + 2v u v z - v z v a + —e 



-3a _ ze -2 a 



(2.8) 



-2a 



+ — e~ 3a )dz 
2(i 



dE L =(2ze~ 2a - ^-e- 3a ^da+(^2vl - 12v u v a + 16^ - 

+ (4zv a - 12zv u Vz^j dv a -\- ^ 12zv a -\- 32zv u -\- 2v z S jdv l 
+ (2v u -v a ^jdv z , 

the secondary constraint $2 can be found as 



iy c dE L = $ 2 =(2v u - v^j (iSzVu - Azv a + 2v z ^ + z^Ae~ 2a - • 



(2.9) 



(2.10) 
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Let us now consider a vector field V which solves the Lagrange field equations 



iruJL = -dE L . 

In light of equations (2.6) and (2.9), equation (2.11) yields by comparison 

T 2 = 2zv a +v z - 2zY a , 
T Vu = \r Va + (2v u - v a )T a +{2v a + ^) (y - v u ) - 1 (V 2 « - je- 3 " 

T Vz = -2zT Va + 4zv a T a - 2v a (2zv a +v z ^+z ^8e" 2a - ^e" 3a ^ , 



(2.11) 



(2.12) 

(2.13) 
(2.14) 

(2.15) 



whereas r a and r„ Q remain arbitrary functions of a, u, z, v a , v u and v z . The dynamics is 
tangent to the constraint manifold if the Lie derivative along V of the secondary constraint 
$2 in equation (2.10) is vanishing, i.e. £r$2 = 0. This means that the contraction i-pd^2 
should vanish, and leads to a restriction on the coefficient r a , which is found to take the 
value T a such that 



(2zv a + v z ^j 


^fe" 3a -2e" 2a + 2(2v u - 


Va) (Sv u - 2v a + ^) 






(8v u -2v a + ^ 





(2.16) 



since the one-form cM>2 is 
'21z 2 



d<& 2 = 



-3a 



8ze 



4v 2 a - 24v u v a + 32v 2 u + 4e" 2a - ^e" 3a ] dz 



+ ^8zv a — 24zv u — 2v z ^jdv ol + (^ — 24zv ol + 64zv u + 4v z ^jdv z 
+ (4v u -2v a S jdv z , 



(2.17) 
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and the various coefficients of r„ Q appearing in i r d&2 add up to zero. The geometrical 
meaning of our calculation is as follows. If Y solves the Lagrange equations (2.11), for any 
vector field Y e ker ujl one finds 

i T+ ~ u L = (i r u^j + (i~ ul^J = -dE L . (2.18) 

In other words, the arbitrariness of r o and Y Va reflects the possibility of adding to any 
solution of equation (2.11) an arbitrary vector field Y G ker ujl- However, if the vector 
field T is also tangent to the constraint manifold, only Y Va remains arbitrary. The vector 
field Yt that solves equation (2.11) and is tangent to the constraint manifold is thus found 
to be 

r - = F 4 + F »li + F 4 +r "-£ + F »-^ + F »-l: ■ (2 - 19) 

where r o has been evaluated as in equation (2.16), and Y u , Y z , Y Vu , Y Vz are values taken 
by the right-hand sides of equations (2.12)-(2.15) at T a = r a . The arbitrariness of Y Va is 
due to the existence of a vertical kernel. 

3. - Generalized method. 

We here study a different and more recent method for the constrained analysis of a 
dynamical system. Interestingly, it does not rely on a Cartan one-form or a pre-symplectic 
two-form as the method studied in section 2, but enables one to derive constraints by 
looking directly at the equations of motion, and can be applied to any system of differential 
equations in implicit form. 

To describe the method, let us assume that a system of implicit dynamical equations 
is given in the particular form 

aij(q,v)v l - fj(q,v) = , (3.1) 
b ij (q,v)(q i -v i ^ =0 . (3.2) 
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One now has a choice of first- or second-order formalism. If first-order theory is used, the 
time-derivatives of the positions q l are not identified with the velocities v % . Denoting by 
ip J a vector such that ip^bij = 0, the most general form of Eq. (3.2) is 

b ij (q,v){q i -v i + ^ =0 . (3.3) 

Thus, in first-order theory, one finds 

9* = ^-^ . (3.4) 

The constraints are found by solving for <tV the equation <^a^ = 0. This leads to a 
compatibility condition of the kind <j>> fj = 0. If this equation is not identically satisfied 
then <j>> fj is the first-generation constraint of the theory. The constraints are preserved by 
requiring that 

In second-order theory one has instead q l = v l , hence t(j l = in (3.4). Note that the 
arguments developed so far do not rely on the existence of a Lagrangian, in agreement 
with what we said. 

However, if a Lagrangian is known for the model under consideration, a^- is the Hessian 

matrix Hij, whereas bij = q^q^ — Qgjg^i ■ Multiplying by a vector A 1 the second-order 
Lagrange equations 

d dL dL 

Jtd^~W = ' ( } 

where 

v\t) = j t q\t) , (3.7) 

one finds 

''.r',1' - , (3-8) 



dv l dqi dq 
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which is a restriction on Cauchy data, i.e. a constraint. Following [17-18], such constraints 
are called first- generation constraints. Clearly, the number of independent eigenvectors cor- 
responding to zero eigenvalues of Hij is N — K, where K is the maximal rank of the Hessian 

matrix. Hence first-generation constraints can be denoted by ip^, m = 1, 2, M < N—K. 
Second-generation constraints (if any) are then obtained by requiring that the evolution 
of the system should be tangent to the first-generation constraint manifold, i.e. 

^=^" = . (3.9) 

Analogously, third-generation constraints are determined by imposing that the evolution 
should be tangent to the manifold defined by first- and second-generation constraints, and 
so on. 

A similar analysis can be repeated for first-order Lagrange equations. This leads to 

A first-generation constraint has been thus obtained within the first-order formalism. All 
further constraints (i.e. second-generation, third-generation and so on) are then found by 
imposing that the evolution of the system should be tangent to the constraint manifold. 

We now apply this new method to our cosmological model. The first-order equations 
of motion are 

2(u-v u y )-(ct-v a > )=0 , (3.11) 
2z (z - v^j - 12z (a - Vo^j + 32z (u - v u ^j = , (3.12) 

4z(a-Vc^j -12z(u-v u ^-(z -v z ^) =0 , (3.13) 

i) z — 6zv a — 6zi) a + 16zv u + 16v u z = , (3-14) 

3 z 2 

-v z + 4zv a + 4zv a - I2zv u - \2zv u + 2ze~ 2a - - — e" 3a = , (3.15) 

4 fi 
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-v a + 2v u 



-2a 



2v 2 a + 12v u v a - 16^ + ± e 
= (^v a - I2v^j (a-Va^j 



-3a 



(3.16) 



Since the Hessian matrix Hij is of rank 2 in our case, its kernel is one-dimensional. The 
vector of the kernel of is found to be 



-Az 



(3.17) 



The equations of motion (3.14)-(3.16) can be seen as a row vector which, multiplied by A 1 , 
yields the first-generation constraint. Such a constraint is thus defined by 



,v d 2 L ■ .-dL d 2 L 

W = A \ . . . q 3 - A 1 — - A 1 —-— 
ov l oq J oq % ov J oq l 



(q J — v jA j 







(3.18) 



Bearing in mind equations (3.11)-(3.13), one finds 



ii? = (2v u -v a )(8zv u -2zv a + v z > ) + ^2e" 2a - ^e" 3 ^ . (3.19) 



The constraint ip l is preserved if = 0, which implies 



(2zv a + v z ) e" 3 - - 2e~ 2a + 2 (2v u - v a ^j (8v u - 2v a + ^) 



a = 



49z- 
4M 



. e -3a _ lQ ze -2a + 4z(2v u - U a ) (8v u - 2v a + ^ 



(3.20) 



By imposing that the evolution of the system should be tangent to the constraint manifold 
we have derived a formula for d, which is the analogue of T a in the Gotay-Nester analysis 
of section 2 (Eq. (2.16)). Thus, by using the generalized method, the constraint analysis 
is completed once we determine d, which remains arbitrary if one looks simply at the 
equations of motion which govern the dynamics. 
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It is now very instructive to perform the constraint analysis within the second-order 
formalism. The second-order Lagrange equations (3.6)-(3.7) read 



v a = a , 



v u = u 



v z = z 



v z - 6v z v a - 6zv a + 16v z v u + 16v u z = , 



—v z + 4v z v a + 4zi) a — 12zv u — 12v z v u + 2ze 



-2a 



4 /i 



,-3a 



o , 



-v a + 2v u - e~ 2a - 2vl + 12v u v a - 16u* + — e" 3a = . 



(3.21) 

(3.22) 
(3.23) 
(3.24) 

(3.25) 
(3.26) 



The kernel of the Hessian matrix is again given by (3.17), and the first-generation constraint 
(3.8) turns out to coincide with the first-generation constraint (3.19). The search for 
second-generation constraints leads to 



a 



49z' 

Til 



e" 3Q - 16ze~ 2a + 4z(2u - aj (8ii - 2a + -^j 



- (2za + zj J^ e ~ 3a - 2e ~ 2a + 2 ( 2 « - «) ( 8 « - 2 " + ^) 



(3.27) 



Requiring that the evolution of the system should be tangent to the constraint manifold, 
and defining 



D = 16ze~ 2a - ^- z z e-* a + 4z{ a - 2u ) ( Sii - 2n -f - 



49 



.2 -3a 



(3.28) 
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one finds 

Da = a 



-3a 



4// 



( - 123z 2 a + 98zz - 48z 2 ?i) 



+ e" 2a (l6za - 16i + 32zw) 

_ V 8 t t -2d+^)(V 2c *--e 
2 V \ /U 

+ 4^2-u - a) I (2a -8u)z-^j 

-3a 



-3a 



: | 8e~ 2a - -e~ 3a 



llz 



2e" 2a + 2 (2w - d) (8u - 2d + 



+ 



4 (2za + zj(a- 2uj (sii - 2d + 

— In, 

(2za + zj ^27zd - Uz^j 



4^ 



z 



+ 2e" 2a f2^d + i 

V / z 



(3.29) 



Note that a substantial difference occurs with respect to a first-order formalism, since 
we find a second-generation constraint which is absent in the first-order case. Hence the 
variable v a is determined, which remains instead arbitrary in the first-order analysis. 

The occurrence of the second-generation constraint, which has no equivalent in Dirac's 
Hamiltonian treatment [2] of section 1, is not a peculiar property of the generalized method 
of this section, but rather of the second-order formalism. A second-order Gotay-Nester 
analysis yields the same result. Indeed, for any point Lagrangian in a second-order theory 
with constraints, the contraction of the energy one-form with a vector field A in the kernel 
of the pre-symplectic two-form yields the constraint 



<t>(q,q) = A s 



d 2 L 



8L 



QqsQqh Qgi 



(3.30) 
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which coincides with the second-order first-generation constraint (3.7)-(3.8). The Lie 
derivative of 4>{q, q) along the second-order vector field 



r = «"^ + r W ' <3 ' 31) 



leads to a further constraint which, in our paper, coincides (by construction) with the 
second-generation constraint (3.27). 

Indeed, our model is not the only case of point Lagrangian where second-order for- 
malism leads to further constraints with respect to first-order formalism. For example, if 
one studies the point Lagrangian (cf. [19]) 

1 ( \ P ( \ 2 

L = -vj + V!q 2 +[l - a)v 2 qi + -(qi - q 2 J , (3.32) 

with a 2 — f3 7^ 0, first-order formalism only leads to the constraint 

<£> 2 = -av 1 +p(q 1 -q 2 S J , (3.33) 

whereas second-order formalism also leads to the further constraint 

<^ 3 =(a 2 -(3)q 2 +(3q 1 -a(3{q 1 -q 2 ^ . (3.34) 

In the latter case, the vector field solving the Lagrange field equations (2.11) takes the 
form 

r=i lr +<i2 4 + (- Q<i2+,3 ( < "-*))|r ■ (3 - 36) 

The additional constraint (3.34) is obtained since velocities have been taken to be time- 
derivatives of positions. 

A further example of the inequivalence of first- and second-order theory is given by 
the following Lagrangian on the tangent bundle of SU{2): 

L = Tr (0-3S- 1 , (3.36) 
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where S G SU(2). The corresponding pre-symplectic 2-form is found to be 

u L = -Ttas^S^dS AS^dS^j , (3.37) 

and the first-order dynamics is given by the vector field X 3 associated with the one- 

it 

parameter subgroup e~° 3 . In second-order formalism, however, there is no dynamics 
compatible with the Lagrangian (3.36). This Lagrangian is a Chern-Simons Term, and 
hence is relevant for modern field theory. 

Interestingly, we may use the constraints (3.19) and (3.27) to derive physical predic- 
tions. In other words, by imposing the first-generation constraint i/j 1 = we may cast the 
second-generation constraint in the form 

= a (^e" 3a - 12^e" 2 ^ + z (de^ - ^e" 3 ^ = . (3.38) 

Equation (3.38) leads to 



a a _ 9z_ 
dz V ^ 



z 



2e a - 



da . (3.39) 



Thus, defining 

x = e" a , (3.40) 

xz = r){Z{x)) , (3.41) 

e(x) = log(x) , (3.42) 

'4- |2 

/(//) = --) - (3.13) 

9 _ M 

differentiation of (3.41) with respect to x and comparison with (3.39) leads to the differ- 
ential equation 

^=v(l + f(v)) • (3-44) 
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The corresponding solution for x may be expressed in the form (see Eq. 1. 126 on page 316 
of [20]) 

By virtue of (3.43), Eq. (3.45) implies, after performing some standard integrals, that 

x 2 z(^xz-2j =A , (3.46) 

where A is an integration constant. Thus, by using the definition (3.40) and the relation 
between z and the scalar curvature: z = ^e 3a R, (3.46) leads to the second-order algebraic 
equation for the scalar curvature 

a 3 R 2 - 8aR - 4A = . (3.47) 

Although it seems impossible to solve by analytic methods the field equations in Hamil- 
tonian (cf. [2]) or Lagrangian form, we have been able to re-express first- and second- 
generation constraints in terms of physically relevant quantities, i.e. the cosmic scale 
factor and the scalar curvature. Interestingly, the roots of (3.47) are given by 



R= ^±^V 4 + Aa • ( 3 - 48 ) 

Thus, if second-order formalism is used, the scalar curvature does not vanish in vacuum in 
the presence of torsion, unless the constant A is set to zero and the negative sign is chosen 
in front of the square root. Moreover, if the constant A is negative, the cosmic scale factor 
is bound to remain less than or equal to 

\A\ 

By contrast, if torsion vanishes, our model corresponds to a Hamiltonian system with 
first-class constraints (cf. [7]). At a Lagrangian level, the first-generation constraint is 
given by El/N, where El is the energy function and N is the lapse function. Thus, such 
a constraint vanishes if and only if the energy function vanishes, and is automatically pre- 
served, since El is constant along solutions of the field equations (hence its Lie derivative 
along the vector field solving the Lagrange equations vanishes) . 
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4. - Concluding remarks. 



The contribution of this paper is a detailed application of the generalized method of 
section 3. Since it relies on an approach [17-18] whose range of applicability is wider than 
any previous (Lagrangian) method, we found it appropriate to focus on this technique 
in our paper. Remarkably, in our specific model, if Lagrange equations are studied in 
second-order formalism, a second-generation constraint is found which is absent in first- 
order formalism. This constraint has been expressed in terms of the physical quantities 
of the problem, as shown in Eq. (3.38). One thus finds the equation (3.47) for the scalar 
curvature, whose roots are given by (3.48). It turns out that the scalar curvature may 
not vanish in vacuum. Interestingly, an upper limit for the cosmic scale factor exists if 
the constant A in (3.48) is negative. Hence second-order theory is found to yield relevant 
information about the early universe in our toy model. 

Note also that the secondary constraint of equation (2.10) is obtained multiplying by 
— \ the secondary constraint of equation (1.7), as one may check by using (2.1) and the 



which expresses the map in local coordinates from the tangent bundle TQ to the cotangent 
bundle T*Q. This implies that, up to an unessential proportionality constant, the Hamil- 
tonian and the Lagrangian method lead to the same secondary constraint. However, the 
careful reader may have noticed that, using tangent-bundle formalism, no primary con- 
straint occurs. Such property is not surprising for the following reasons. If det | d ®ig vj |= 0, 
the map FL defined in (4.1) maps the tangent bundle TQ to the primary-constraint sub- 
manifold E of the phase space T*Q. Thus, (FL)* maps (by definition of pull-back) forms 
on S to forms on TQ, and the forms ul, alE^ appearing in the Lagrange equation (2.11) 
correspond to geometrical objects already living on E C T*Q. Interestingly, the first set of 
constraints one actually evaluates in the Lagrangian formalism are the equivalent of sec- 
ondary constraints of Dirac's theory. This property appears very important for first-class 
theories such as Maxwell's electrodynamics and GR: the Lagrangian method is a powerful 



map 




(4.1) 
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way to evaluate directly the constraints which govern the theory, i.e. first-class secondary 
constraints. Although the question remains open [11-13,21], Dirac's consideration of first- 
class secondary constraints appears more natural within the Lagrangian approach. 

A naturally occurring question is whether a Lagrangian multisymplectic analysis [22] 
of first-class theories such as Einstein's GR offers any advantage with respect to the ADM 
Hamiltonian method [13,23-24] which relies instead on the 3+1 split of the space-time 
geometry, especially in light of the longstanding unsolved problems of the quantization 
program [13]. Work is now in progress by the authors on some of these issues, and we hope 
that the geometrical reformulation of certain properties of classical field theory presented 
in this paper will serve as a first step towards this more ambitious goal. 

* * * 
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